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ABSTRACT  ^ 

> 

~~"*bynam1c  extension  of  Slh’s  fracture  criterion  based  on  strain  energy 
density  factor*  ^(dW/dV)*  Is  used  to  analyze  dynamic  crack  propagation  and 
branching.  Influence  of  the  nonsingular  components*  which  are  known  as  the 
higher  order  terms  (HOT)  in  the  crack  tip  stress  field*  on  the  strain  energy 
density  distribution  at  a  critical  dlstance^sur rounding  the  crack  tip  moving  at 
constant  crack  velocity  Is  examined.  This  r£(dW/dV)  fracture  criterion  Is  then 
used  to  analyze  available  dynamic  photoelastic  results  of  crack  branching  and 
of  engineering  material sw 
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INTRODUCTION 


The  mechanics  associated  with  the  curved  crack  paths  or  bifurcation  events 
received  continuous  attention  In  the  literature.  Comprehensive  reviews  on 
crack  curving  and  branching  were  published  by  the  authors  Cl]  and  Dally  C23. 
Literature  on  crack  branching  during  rapid  fracture  discusses  Its  correlation 
with  either  the  dynamic  crack  tip  stress  field*  Initiation  of  secondary  cracks* 
or  with  stress  wave  loading.  These  studies  deal  primarily  with  the  stress 
Intensity  factor*  the  strain  energy  release  rate  at  the  Instant  of  branching 
and  the  branching  crack  velocity.  The  crack  branching  angle  as  well  as  the 
crack  curving  angle  can  also  be  determined  by  a  crack  directional  stability 
criteria  under  combined  tension  and  shear  loadings*  l.e.*  for  a  mixed  mode 
stress  field.  This  crack  stability  criterion  Is  governed  by  either  the  maximum 
circumferential  stress  at  the  crack  tip  C3]  or  at  a  specified  distance  from  the 
crack  tip  [4-6]*  the  maximum  energy  release  rate  [7*8]*  and  the  minimum  strain 
energy  density  factor  surrounding  the  crack  tip  [9*10]  or  at  a  specified 
distance  from  the  crack  tip  [11].  Slh  [12]  extended  the  minimum  strain  energy 
density  factor  criteria  to  rapidly  propagating  cracks  for  predicting  the  crack 
branching  angle  with  a  singular  stress  field. 

Rossmanlth  [13]  and  the  authors  [14-16]  have  shown  that  the  nonsingular 
component  In  the  stress  acting  parallel  to  the  crack  path  has  considerable 
Influence  on  the  directional  stability  of  a  moving  crack  and  Incorporated  this 
component  for  estimating  the  curving  and  bifurcation  angles.  The  additional 
higher  order  terms  In  the  crack  tip  stresses*  which  are  used  to  account  for 
the  spread  of  advanced  cracks*  are  the  cause  of  Incipient  branching  In  brittle 
materials*  and  enlarge  the  fracture  process  zone  which  results  In  roughening  of 
the  fracture  surface.  The  presence  of  unsuccessful  multiple  cracks  or  crack 
Interaction  with  the  stress  waves  [17-19]  under  mode  I  loading  are  often 
associated  with  the  slight  unsymmetry  In  the  fringe  patterns  of  the  dynamic 
photoelastic  experiments  of  rapidly  propagating  cracks.  This  unsymmetry  In  the 
fringe  pattern  was  modeled  by  Rossmanlth  [17]  and  the  authors  [19]  by  adding 
the  higher  order  terms  of  mode  I  and  mode  II  stress  components  to  the  mode  I 
singular  crack  tip  stress  field.  This  inevitable  Involvement  of  the  higher 
order  terms  forms  the  basis  of  the  dynamic  fracture  analysis  In  this  paper  and 
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the  prediction  of  crack-growth  direction  of  a  constant  crack  velocity  using  the 
strain  energy  density  factor*  rc(dW/dV)»  criterion. 

DYNAMIC  MIXED  MODE  CRACK  TIP  STRESS  FIELDS 


The  crack  tip  dynamic  state  of  stresses  under  mixed  mode  conditions  was 
given  by  Nlshloka  and  Atlurl  [20]  In  terms  of  the  local  rectangular  (x»y)  and 
polar  (r*6)  coordinates.  The  three  rectangular  stress  components  under  mode  I 
and  mode  II  conditions  are  given  as: 
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r  and  0  are  the  polar  coordinates  with  the  origin  at  the  moving  crack  tip#  C» 
C^#  and  C2  are  the  crack  velocity#  dllatatlonal  and  dlstortlonal  wave 
velocities#  respectively. 

The  general  solution  expressed  In  Equation  (1)  yields  the  singular 
stresses  when  n  »  1,  l.e. #  A^  »  Kj  and  *  Kjj#  which  are  stress  Intensity 
factors  of  mode  I  and  mode  II#  respectively.  The  constant  stress#  a Is 
related  to  the  higher  order  term  (HOT)  for  n  =  2. 

6Bt(c)(S?-S?) 


DYNAMIC  STRAIN  ENERGY  DENSITY  FACTOR  CRITERION 


The  strain  energy  density  factor#  rc  (dW/dV)»  criterion  developed  by  SI h 

[9]  Is  based  on  the  local  density  of  the  strain  energy  field  at  a  core  radius 

rc  from  the  crack  tip.  It  assumes  that  fracture  Initiates  when  rc(dW/dV) 

reaches  a  critical  value  of  [r_(dW/dV]„. 

c  c 
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The  elastic  strain  energy  density*  W,  per  unit  volume*  V*  of  the  material 


Is: 


■jjj  »  C(0  2  +  O  2  +  o  2)  - 

dV  2E  xx  yy  zz 


T—  (O  +o  +  o  )2  +  20  2]  (3) 

l+v  xx  yy  zz  xy 


for  plane  strain  condition: 
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where  E  and  v  are  the  modulus  of  elasticity  and  Poisson's  ratio*  respectively. 


As  discussed  earlier*  for  n  *  2  the  first  nonsingular  term  exists  only  In 
oxx*  stress  component  of  mode  I  loading  and  Is  zero  In  all  other  components  of 
mode  I  and  mode  II  crack  tip  stress  fields.  These  and  other  higher  order  terms 
for  n  >  2,  on  the  strain  energy  density  factor  surrounding  the  crack  tip  are 
believed  to  Influence  the  spreading  of  advanced  cracks  which  turn  leads  to 
crack  curving  and  branching.  By  denoting  the  energy  density  function  (dW/dV) 
In  the  presence  of  higher  order  terms  as  (dW/dV>n*  the  Intensity  of  the  strain 
energy  density  factor  for  the  state  of  plane  strain  at  a  critical  distance*  rc» 
can  be  written  for  as:  rc(dW/dV>n.  The  subscript  Index,  n.  Indicates  the 
number  of  the  terms  considered  In  determining  the  strain  energy  density  factor. 
By  using  Equations  (1)  through  (5)  and  evaluating  for  nsl*  the  strain  energy 
density  factor*  S*  became: 


s 


Call  AIl' 


+  2a12  An 


AIIX  +  a22 


(6) 


The 


The  coefficients  a^  depend  upon  C»  C^»  C2  and  Poisson’s  ratio*  v  . 
fracture  criterion  based  on  the  minimum  strain  energy  density  factor  notes  that 
the  crack  will  extend  to  the  location  of  the  minimum  strain  energy  density 
factor#  Crc(dW/dV)]mjn.  In  the  presence  of  higher  order  terms#  strain  energy 
density  factor*  criterion  can  be  written  as: 
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The  half  crack  branching  angle#  0C#  at  which  CrcdW/dV]n  possesses 
stationary  values*  can  be  obtained  from  Equation  (7).  As  a  check  on  the 
accuracy  of  Equation  (5)  and  (7)*  the  dynamic  r  (dW/dV)  should  coincide  with 

C 

Its  static  value  In  Reference  [123#  when  the  crack  velocity  c — >  0  for  n*l. 


NUMERICAL  RESULTS 


In  order  to  study  the  variation  of  the  strain  energy  density  factor  with 
angular  orientation#  the  theoretical  (dN/dV)  values  were  computed  for  each  0.25 
degrees#  using  a  Poisson's  ratio#  v  *  0.36#  Youngs  Modulus*  E  *  3.74  GPa# 
critical  radius  rc  »  1.3mm#  Aj^  ■  Kj  *  1  MPa<4i»  and  Kjj/Kj  *  0  value.  Static 
strain  energy  density  factor  was  approximated  by  setting  C  =  0.01C2  In 
Equations  (1)  and  (7)  where  the  results  agreed  well  with  that  of  References 
[12*13]. 


Crack  velocities  of  C/C2  *  0.01  to  0.8  were  chosen  to  generate  the  dynamic 
strain  energy  density  factor  In  the  vicinity  of  the  crack  tip.  The 
representative  static  and  dynamic  strain  energy  density  factor  distributions  In 
the  vicinity  of  the  crack  tip  at  a  distance  r  *  1.3  mm  are  shown  In  Figure  1 

C 

for  A^<0»  AI2  «  0#  and  Aj2>0  under  mode  I  loading. 

Figure  2  shows  the  angular  variation  of  static  and  dynamic  (r  dW/dV)  with 
the  third  higher  order  term  Aj^  which  Is  compressive  In  nature*  In  addition  to 
remote  stress  under  pure  mode  I  loading  (l.e.#  Kjj  *  0).  Figure  2(a)  shows  the 
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strain  energy  density  factor  distribution  for  Aj2  ■  -0.25  and  Ajj  *  -0.3  for 
v=  0.36  and  Figure  2(b)  represents  the  distribution  of  strain  energy  density 
factor  for  A^  *  0.25  and  A^  *  -0.3.  The  minimum  value  of  strain  energy 
density  factor  Is  almost  constant  within  +10°  for  A^  <  0  at  lower  crack 
velocities.  However*  the  minimum  strain  energy  density  factor  for  A^  >  0  at 
lower  crack  velocities  yields  minimum  r  (d#/dV)  at  greater  than  ±  30°  when 
crack  velocity  C/C2  >  0.5.  This  Is  not  surprising  since  the  compressive  remote 
stress  stabilizes  the  crack  path  and  the  tensile  stress  enhances  the 
directional  Instability  of  a  propagating  crack. 

Figure  3  shows  the  angular  distribution  of  dllatatlonal  and  dlstortlonal 
strain  energy  density  factors  for  the  strain  energy  density  factors  In  Figure  1 
for  Aj^=  1.0.  It  can  be  observed  that  the  energy  Is  constant  between  ±13 
degrees  for  lower  crack  velocities  and  the  maximum  magnitudes  of  the 
dlstortlonal  energy  and  corresponding  dllatatlonal  energy  are  minimum  at  C/C2  i 
0.5.  The  minimum  or  maximum  spread  of  dllatlonal  and  dlstortlonal  energy  Is 
constant  within  the  +13  degrees  angular  region#  respectively.  A  similar 
observation  was  made  when  we  added  A^  term  to  the  stress  field. 

Figure  4  shows  the  branching  angle  versus  crack  velocity  corresponding  to 
a  singular  stress  field  with  varying  Poisson's  ratio.  The  minimum  strain 
energy  density  occured  at  angle  6*0  for  crack  velocities  C/Cj^  £  0.46;  It 
becomes  non-zero  for  C/C ^  >  0.46.  Note  that  the  branching  angle  Increased  with 
Increasing  crack  velocities.  However#  at  higher  velocities#  the  difference  In 
branching  angle  Is  very  small  regardless  of  the  difference  In  magnitudes  of 
Poisson's  ratio. 

Table  1  shows  the  predicted  crack  branching  angles  using  only  the  singular 
stress  field  with  critical  crack  velocities  of  various  materials  represented  by 
Poisson's  ratio.  These  results  suggest  that  the  branching  angle  Is  very 
sensitive  to  the  Poisson's  ratio#  v.  Note  that  the  branching  angle  of  ill  to 
+23  degrees  occurred  In  these  materials  within  the  crack  velocity  range  of 
C/C2  ■  0.47  to  0.64#  which  Is  higher  than  the  experimentally  observed  crack 
velocities  C13  Irrespective  of  the  material  used.  The  discrepancy  between  the 
experimental  observation  and  theory  can  be  partly  attributed  to  the  finite 
geometry  and  boundary  effects  which  Influences  all  singular  and  higher  order 
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terms  In  the  crack  tip  stress  field.  In  the  following*  some  representative 
values  of  the  higher  order  terms  (HOT)  are  used  to  show  the  difference  In  the 
predicted  branching  angles  with  HOT  and  without  HOT  In  the  stress  field. 

Figure  5  shows  the  branching  angle  versus  the  crack  velocities  for  v® 
0.33*  In  the  presence  of  higher  order  terms  under  the  pure  mode  I  loading 
condition.  The  addition  of  positive  higher  order  terms  enhanced  the  branching 
angle  while  negative  higher  order  terms  suppressed  branching.  The  crack  Is 
directionally  stable  even  at  higher  crack  velocities  of  C/C2  =  0.6  In  the 
presence  of  negative  HOT  terms.  Also*  at  higher  crack  velocities*  the 
resulting  branching  angle  Is  not  Influenced  by  the  higher  order  terms,  with 
differences  of  nearly  8  degrees. 

Quite  often  the  dynamic  photoelastic  fringe  patterns  associated  with  a 
propagating  curved  crack  exhibits  a  slight  unsymmetry  In  Isochromatics  under 
mode  I  loading.  This  unsymmetry  In  Isochromatics  Is  modeled  by  adding  the  mode 
II  higher  order  terms*  without  the  model  II  singular  term.  Since  the  applied 
remote  load  Is  only  mode  I  loading*  It  Is  appropriate  to  disregard  the  singular 
term  and  add  only  AjJ3  term  as  was  demonstrated  previously  In  References 
[18*19]. 

Figure  6(a)  shows  the  angular  distribution  of  the  strain  energy  density 
factor  for  crack  velocities  of  C/C2  *  0.01  to  0.8  and  AI2  *  -0.25*  Aj3  *  0.3* 
and  AII3  *  -0.3  for  v  =  0.36.  The  rc(dW/dV)3  distribution  Is  slightly 
unsymmetrlc  at  lower  crack  velocity  and  Increased  at  higher  crack  velocities. 
For  a  crack  velocity  of  0.4  <  C/C2  >,  0.7,  the  minimum  strain  energy  density 
factor  Is  located  at  about  20°  but  at  C/C2  *  0.8*  the  minimum  strain  energy 
density  factor  occurs  as symmetrically  at  about  55  and  60°.  Reversing  the  sign 
of  AII3  from  negative  to  positive  In  Figure  6(b)  simply  reverses  the  peaks  In 
the  strain  energy  density  factor  distribution.  At  lower  crack  velocities  the 
positive  higher  order  terms  of  mode  II  stress  field  enhances  crack  kinking  and 
unsymmetrlc  crack  branching  could  occur  at  higher  crack  velocities.  Therefore, 
for  AII3  <  0*  a  positive  kinking  angle  Is  obtained  while  for  Ajj3  >  0*  a 
negative  kinking  angle  results.  This  effect  Is  similar  to  the  presence  of 
positive  and  negative  mode  II  stress  Intensity  factors,  respectively.  The 
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series  of  experimentation  shows  that  branching  could  occur  even  In  the  presence 
of  Kjj  at  higher  crack  velocities. 

Qualitative  results  of  crack  branching  angles  are  presented  In  the 
following  for  Poisson's  ratio  of  v  =  0.25#  0.29#  and  0.33  for  =  1.0#  Aj2  * 
-0.25#  Ajg  =  0.3  and  higher  order  terms  Involving  mode  II  nonsingular  terms  of 

AII1  *  AII2  =  an<J  AII3  =  °*3*  Tabl®  2  shows  the  effect  of  HOT  on  crack 
branching  angles  for  varying  crack  velocities.  At  lower  crack  velocities  of 
C/C2  <  0.6#  the  mode  II  higher  order  terms  contributed  to  crack  kinking  and 
branched  asummetrlcal 1y  at  higher  velocities.  This  observation  Is  consistent 
with  the  present  results  In  that  the  HOT  effect  Is  negligible  at  higher  crack 
velocities.  An  enhanced  branching  angle  also  Is  clearly  observed  at  higher 
velocities  Irrespective  of  the  number  of  higher  order  terms  present.  Positive 
higher  order  terms  enhanced  magnitudes  of  the  branching  angle  at  all  crack 
velocities.  However,  the  crack  branching  angle  was  always  0  degree  when  C/C 2  < 
0.4  for  all  compressive  or  negative  higher  order  terms. 

Representative  experimental  crack  branching  data  associated  with  running 
cracks  In  Homal lte-100#  polycarbonate#  steel#  aluminum  [21]  and  glass  were 
evaluated  and  presented  In  Table  3.  There  Is  good  agreement  between  measured 
and  predicted  angles  using  higher  order  terms#  whereas  significant  differences 
are  observed  when  HOT  terms  were  not  used. 

CONCLUSIONS 

(1)  The  minimum  value  of  rc(dW/dV)  strongly  depends  on  the  first  nonsingular 
stress#  0Qx#  for  a  given  Poisson's  ratio  under  any  mode  of  deformation. 
It, also  depends  strongly  on  the  crack  velocity  and  varies  systematically 
with  the  nonsingular  stresses. 

(2)  For  small  values  of  C/C2  <  0.4#  the  fracture  angle  ec#  Is  found  to  be  0 
for  all  negative  higher  order  terms  of  mode  I  loading  and  ec  ^  0  for  all 
positive  values  of  HOT#  Irrespective  of  crack  velocity.  Negative  higher 
order  terms  suppresses  the  crack  extension  angle  0  whereas  positive  HOT 

G 

always  enhances  the  0c#  regardless  of  the  Poisson's  ratio. 
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(3)  Depending  on  the  test  conditions  and  the  shape  of  the  fractured  specimens* 
theoretically  predicted  angle  could  deviate  from  experimentally  measured 
fracture  angles  when  the  second  order  term  Is  neglected. 
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TABLE  2 


EFFECT  OF  MODE  II  HIGHER  ORDER  TERMS  ON  CRACK 
EXTENSION  ANGLE 


An  *  1.0#  AI2  =  0.25#  AI3  *  0.3  AND  An3  *  0.3« 


Pol sson ' s 

Crack  Velocity 

Crack  Extension 

Remarks 

Rati o  v 

c/c2 

Angle  ©c 

0.50 

2.50 

Crack  kinking 

0.54 

11.75 

Crack  kinking 

0.58 

33.25 

Crack  kinking 

0.25 

0.62 

±  44.25 

Branching 

0.66 

±  51.00 

Branching 

0.72 

+  57.75 

Branching 

0.56 

5.25 

Crack  kinking 

0.60 

27.00 

Crack  kinking 

0.29 

0.64 

+  40.75 

Branching 

0.70 

+  52.00 

Branching 

0.58 

3.0 

Crack  kinking 

0.33 

0.64 

30.25 

Crack  kinking 

0.70 

+  47.00 

Branching 

Crack  extension 

angles  are  negative  for 

*113  < 

TABLE  3 


MEASURED  AND  CALCULATED  VALUES  OF  BRANCH  ANGLE 


Material 

Crack 

r. 

Measured 

Predicted  Remarks 

Velocity 

mm 

Branch  Angle 

Angle  2  0 

c/c2 

2  % 

Without  Hot  WTth  Hot 

Steel 

0.44 

1.0 

66 

0 

64 

Pressurized 

pipe 

Aluminum 

0.14 

1.3 

88 

0 

84 

Pressurized 

pipe 

H0MALITE-100 

0.10 

1.3 

52 

0 

52 

WL-RDCB 

specimen 

0.46 

1.3 

50 

0 

47 

B1  axially 

loaded 

Centrally 

notched 

specimen 

0.38 

1.3 

28 

0 

0 

SEN 

specimen 

Pol ycarbonate 

0.45 

0.75 

34 

0 

0 

SEN 

specimen 

DISTORTIONAL  STRAIN  ENERGY  DILATATIONAL  STRAIN  ENERGY  DENSITY  FACTOR 

DENSITY  FACTOR 


Au  =  1.0 


C/C2  -  0.8 


-29.0  20.0 


100.0  U0.0  10 


C/0 ;  =  0.8 


100.0  -11*0.0  -100.0  -00.0  -20,0  20.0  00,0  l#0>#  #  u 

Al'lGIf,  8  (b) 

FIGURE  3.  ANGULAR  DISTRIBUTION  DILATIONAL  AND  DISTORTIONAL  STRAIN 
ENERGY  DENSITY  FACTORS  UNDER  SINGULAR  STRESS  FIELD. 


FIGURE  EFFECT  OF  POISSON'S  RATIO,  ,  ON  CRACK  BRANCHING 

ANGLE  UNDER  NODE  I  SINGULAR  STRESS  FIELD  CONDITION. 
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